Enhancement of corona discharge induced wind generation with carbon nanotube and titanium dioxide decoration
Ye Jianchun1, Li Jun2, Chen Xiaohong1, Huang Sumei1, Ou-Yang Wei1, ‡
Engineering Research Center for Nanophotonics and Advanced Instrument (Ministry of Education), School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
Department of Electronic Science and Technology, Tongji University, Shanghai 201804, China

 

† Corresponding author. E-mail: ouyangwei@phy.ecnu.edu.cn

Abstract
Abstract

Dip-coated double-wall carbon nanotubes (DWCNTs) and titanium dioxide (TiO2) sol have been prepared and smeared onto the tip of a conductive iron needle which serves as the corona discharge anode in a needle–cylinder corona system. Compared with the discharge electrode of a CNT-coated needle tip, great advancements have been achieved with the TiO2/CNT-coated electrode, including higher discharge current, ionic wind velocity, and energy conversion efficiency, together with lower corona onset voltage and power consumption. Several parameters related to the discharge have been phenomenologically and mathematically studied for comparison. Thanks to the morphology reorientation of the CNT layer and the anti-oxidation of TiO2, better performance of corona discharge induced wind generation of the TiO2/CNT-coated electrode system has been achieved. This novel decoration may provide better thoughts about the corona discharge application and wind generation.

1. Introduction

Corona discharge-induced wind (CDIW), simply called ionic wind or electric wind, is an airflow induced by corona discharge which is generated from a discharge electrode with a sharp tip (such as a needle, a wire, or a blade) under high voltage. The ionized molecules around the tip accelerate under the high electric field and collide with the surrounding neutral molecules in air, forming an airstream between the electrodes.[1] Compared with a traditional mechanical fan composed of a motor, a speed-regulator, and flabella, CDIW has at least three advantages. (i) It is operated with an extremely low electrical current ( ), indicating that CDIW is extremely energy-saving and its energy consumption is two orders of magnitude lower than that of the traditional mechanical fans.[2] (ii) CDIW instruments can be fabricated in small size without any moving parts. (iii) CDIW can be easily and quickly generated in atmospheric air, making it suitable for fast-response applications. Hence, CDIW has recently received lots of attention in the fields of heat transfer for cooling devices,[3,4] drying and thawing food,[5,6] and other applications such as electrostatic precipitators,[7] plasma actuating systems,[8] air purification,[9] and electrohydrodynamic thrusters.[10,11] In these applications, the strength of CDIW is a prerequisite and thus strategies to improve CDIW are of particular importance.

Most CDIW instruments consist of a needle-to-collecting electrode system, where the needle with a smaller diameter can generate a large local field intensity when other parameters are fixed. Numerous experimental and analytical investigations have been conducted to study the corona discharge behavior with respect to a series of parameters, such as various electrode configurations,[12,13] wind velocity profile,[1416] voltage polarity,[16,17] air moisture,[9,14,17] as well as electrode gap.[12,15,18] Meanwhile, nanoscale materials with inherently high aspect ratios, small tip radii, and sharp edges, such as carbon nanotubes (CNTs),[19,20] graphene sheets,[21] ZnO nanowires,[22] and other nanostructures,[2325] have been investigated for the development of stable atmospheric corona discharges. Among the nano-materials, a thin film of CNTs has been demonstrated as an ideal cathode candidate for electron emission in field emission devices due to the high aspect ratio, outstanding electrical conductivity, excellent chemical stability, and superior mechanical strength and durability.[20,26,27]

However, few studies have previously been carried out to investigate the CDIW of a sharp discharge electrode combined with nanoscale materials from the viewpoint of improving the strength of CDIW. Recently, we invented a simple method to decorate the tip of a discharge needle electrode with a dip-coated CNT thin film and found that the strength of CDIW was greatly improved.[2] Meanwhile, our recent studies on the effect of a thin film of titanium dioxide (TiO2) nanoparticles coated on the CNT cathode of the field emission devices have presented outstanding properties due to the enhanced edge effect arisen from the unique TiO2/CNT micro films.[28,29] Inspired by the synergistic effect of the TiO2/CNT composite, in this article, we utilized dip-coated double-wall carbon nanotube (DWCNT) slurry and prepared TiO2 sol to decorate a conductive iron needle tip as the discharge electrode so as to acquire better efficiency of ionic wind.

2. Experiments

An ordinary iron needle was used as the referential discharge electrode with a diameter of 2.54 mm measured by the thickness gauge instrument. The bare needle was cleaned following a common cleaning process for silicon wafers[30] by acetone, isopropanol, and deionized water sequentially under an ultrasonic environment. For comparison, two kinds of other discharge electrodes were made based on the bare electrode. One was prepared by coating a thin layer of CNTs onto the tip of the needle, and the other was fabricated by further coating a thin layer of TiO2 on the surface of the CNT thin film. For preparation of the CNT-coated electrode, DWCNTs (purity , diameter , length , Shenzhen Nanotech Port Co., Ltd, China) were used and the details of the coating process were described in our previous article.[2] For the fabrication of the TiO2/CNT-coated electrode, TiO2 sol was first synthesized through the hydrolysis and condensation of tetrabutyl titanate (Ti(OC4H9)4) following our previous article and then dip-coated onto the CNT-coated discharge electrode following by a drying process of 200 °C for 10 min and a subsequent annealing process of 350 °C for 30 min.[28]

The experimental appliance is sketched in Fig. 1. The iron needle served as the discharge anode while the thin copper cylinder with an inner diameter of 19.46 mm was the cathode. The needle electrode was mounted at the axis of the cylinder electrode to achieve optimum performance and the distance between the electrodes (also called the electrode separation or electrode gap) was fixed at 5.5 mm. The power was provided by a high-voltage direct current power supply which could generate and regulate the output voltage from 0 to 30 kV with the resolution of 0.01 kV. The electric current was measured using a multimeter with the accuracy of . The needle was connected to the negative pole of the high voltage power supply while the copper cylinder was safely grounded. The whole corona wind generating device was positioned in a polymethyl methacrylate (PMMA) tube with the inner diameter of 19.50 mm. The hot wire anemometer (AM-4204HA, Lutro Electronic Enterprise Co. Ltd, China) with the range of and resolution of was attached to the outlet of our apparatus to examine the practical velocity of the ionic wind. The corona discharge test was conducted at ambient air. Note that the measured wind velocity at the outlet is larger than the average speed of the airflow in the PMMA tube and it can be amplified by assembling a number of discharge electrodes in parallel.[1]

Fig. 1. Sketch of experimental appliance for CDIW generation. The diameter of the copper cylinder electrode is 19.46 mm, the inner diameter of the PMMA tube is 19.50 mm, and the distance between the electrodes is 5.5 mm.
3. Results and discussion

The relationship between the corona current (I) and the applied voltage (V) for various experimental electrodes is plotted in Fig. 2. All three curves in Fig. 2(a) tend to become steeper as the voltage increases. This means that the electrical resistance of the gas channel between the two electrodes tends to lessen when a higher voltage is applied especially as the discharge devices are close to spark over. In addition, once the applied voltage reaches the corona onset voltage (V0, defined at the corona current of ), the electrical current flows from the needle to the grounded electrode because of the ionization of air near the discharge needle tips. The onset voltage is the smallest for the TiO2/CNT-coated electrode (5.9 kV) whilst the largest for the bare electrode (6.9 kV). As a result, when an identical corona voltage is applied, the TiO2/CNT-coated device can acquire the highest discharge current. The diminutive discharge onset voltage is necessary and desirable for industrial applications of CDIW devices as low V0 probably diminishes the power consumption thus increases the energy conversion efficiency.[13] The helpful ways of minimizing V0 can be adopted by reducing the discharge size as well as the inter-electrode spacing. However, the lowered electrode separation may cause reduced ionic wind velocity while the identical input power is exerted.[31,32] As a consequence, external decoration onto the tips/edges of the discharge electrodes is of great significance to obtain minor onset voltage. Since the air between the two electrodes would almost not be heated during the discharge process, the temperature of air can be basically unchanged during the discharge process. As a result, the descended V0 for the TiO2/CNT-coated device is unrelated to the temperature.[33]

Fig. 2. Current–voltage characteristics of the as-prepared electrodes: (a) IV curves, (b) I/VV curves (the red lines are linear fittings for the experimental data).

Figure 2(b) depicts the ordinary conductivity with respect to the applied voltage (I/VV curves) of the as-prepared devices. The almost linear curves manifest that the discharge performance follows the Townsend’s relation for all of our devices[34] where I is the corona discharge current (in A), V is the applied voltage (in V), V0 is the corona onset voltage (in V), and B is the dimensional constant depending on the mobility of charge carriers in the drift region and the geometrical factors of the devices such as the inter-electrode distance, the radius of the needle electrode, etc. The fitting slopes B in Fig. 2(b) are listed in Table 1. Herein, we have recognized that there is another definition form of conductivity, i.e., the differential conductivity, which has been intensively studied in graphene related electronic devices recently.[3537] Therefore, we plot the differential conductivity with respect to the applied voltage for guiding the eyes (see Fig. S 1 in the supplementary material) and the results confirm that the discharge process follows the Townsend’s relation.

Table 1.

Calculated parameters of various electrodes.

.

The slope B in Fig. 2(b) can be elaborately defined to search the ionic mobility (μ) of corona discharge[2] where ε is the absolute dielectric permittivity of the gas medium in , which approximates to the dielectric constant ε0 of vacuum , μ is the carrier mobility of the ionized gas in , and r and R are the radii of the needle and cylindrical electrodes in m, respectively. The calculated ionic mobilities in our experimental devices correspond to the ionic mobility of common discharged air,[1] and are comparable with the state-of-the-art organic and inorganic semiconductors.[38,39] What seems interesting is that the CNT-coated electrode manifests the highest ionic mobility of while the bare electrode reveals the minimum ionic mobility of only . Therefore, DWCNTs and the further coated TiO2 thin layer do have an effect on accelerating the migration of ions.

In the gas phase, the ionic mobility can be defined as where e is the elementary charge in C, me is the electron mass in kg, and vm is the collision frequency for electron momentum transfer in s−1. From Eq. (3), vm can be obtained for various experimental electrodes and listed in Table 1 for comparison. Noted that the calculated collision frequency is three orders of magnitude higher than the reported results because the international system of units is used in this study.[2,40]

The electric wind velocity (v) as a function of the applied voltage (V) is shown in Fig. 3. The wind velocity increases more sharply in the lower voltage region (below 8 kV) and then ascends more slowly in the higher voltage region. The TiO2/CNT-coated electrode can achieve a larger velocity than the other two counterparts when the identical voltage is applied and manifest the maximal peak velocity of . The connection of the current with the ionic wind velocity is given by[34] where g1 is an appropriate function of the geometry and inversely proportional to the square root of the aerodynamic loss coefficient K, and ρg is the density of air (in . The Reynolds number (Re) for Newtonian fluid (e.g., air flow in this work) flowing in a tube can be defined as where D is the diameter of the tube in m and η is the dynamic viscosity of the ionic air in . As we know, the laminar flow is generally determined with while the turbulent flow is defined when for the fluid in a tube. The Reynolds number between 2000 and 4000 is the transition flow from laminar to turbulent flow.[41] Supposing the ionic wind velocity is (which is comparable to our maximal wind velocity available), the simply calculated Reynolds number is 2358. Considering the wind velocity in the tube is much smaller than that recorded at the outlet, the real Reynolds number shall be much lower than the critical value of 2000 and thus the air flow in the CDIW tube is a kind of laminar flow. Furthermore, the air velocity can be even smaller at the lower voltage region which symbolizes that laminar flow plays a leading role for CDIW in a tube. Note that the larger Reynolds number obtained from the TiO2/CNT-coated electrode contributes to the lower aerodynamic loss coefficient K.[2]

Fig. 3. Ionic wind velocity v as a function of discharge voltage V for various electrodes (the red lines are linear fittings for the data at high voltage regions).

Using the following two parameters: and correlating Eq. (2) with Eq. (4), we can obtain

The vV curves (as shown in Fig. 3) in the high-voltage domain are approximately linear, which accord to the deduced equations above. The fitting slopes, i.e., the slope g3 in Eq. (9), of our experimental data for the three electrodes are listed in Table 1 as well and the slope g3 reduces successively for bare, TiO2/CNT-coated, and CNT-coated electrodes. There are several merits of CNTs such as high electrical conductance, glorious mechanical properties, as well as unique aspect ratios.[27,4244] It is convenient to form an extremely high local electric field at the vicinity of the tips. The annealing procedure of the TiO2 layer can produce plenty of teeny modules on the TiO2/CNT nano-composite thin membrane. CNTs become more prominent and upright under the impact of these modules, leading to the diminishment of the shielding effect as well as the enhancement of the fringe effect.[28] The local electric field strength adjacent to the tip of the TiO2/CNT-coated discharge electrode would be even huge contributing to the declined onset voltage and incremental wind velocity. Moreover, the tendency of the vV curves for our electrodes approaches to the computational fluid dynamics results and thus the energy loss is less.[45]

The power consumption as well as the operation efficiency may be of great concern for electrical devices. It is urgent to search brilliant low-powered CDIW cooling apparatus in order to replace conditional cooling fans. The kinetic output power P0 (in W) as a function of CDIW velocity v is given as[34] where S is the cross-sectional area of the PMMA tube at the outlet (in m2) of the PMMA tube for our experimental setup. The calculated kinetic output power P0 using the above equation depending on CDIW velocity v is given in Fig. 4(a). Among our fabricated electrodes, the TiO2/CNT-coated one gives the optimal output power. The electrokinetic energy conversion efficiency (ECE) ηe is defined as[34] where Pi is the input power in W.

Fig. 4. (a) Corona discharge output power (P0) and (b) efficiency (ηe) as a function of voltage (V) for various electrodes.

The electrokinetic energy conversion efficiency is commonly knockdown for CDIW devices ( ), impeding the practical application severely.[31,34,46] However, there is scarcely any effective techniques to improve ECE. Higher ECE can be achieved by decreasing the input power (i.e., the power dissipation) and/or increasing the output power. Therefore, the way to attain high output work is quite imperative and the novel electrode decoration in this article can surprisingly give rise to better ECE. As seen in Fig. 4(b), the TiO2/CNT-coated electrode exhibits the best ECE. That is because the ionic wind velocity v of the TiO2/CNT-coated electrode is the largest among the three electrodes when the applied voltage is the same (manifested in Fig. 3). Consequently, the corona discharge performance of the TiO2/CNT-coated electrode is optimal, which makes it a promising candidate for CDIW generators. The decline of ηe at the applied voltage over 10 kV is probably caused by the sharply decreased electrical resistance of the air channel between the two electrodes as the discharge apparatus nearly breaks down which is in accordance with the discussion of Fig. 2(a).

Furthermore, DWCNTs have poor adhesion to the conductive substrate, causing an unwanted contact resistance, which in turn creates the heating effect.[4749] The poor adhesion may be a barrier for electrical carriers to transport during the discharge process. The created heating effect would ascend the electrical resistance of the discharge needle and then the corona current falls off. DWCNTs can be easily degraded as a result of oxidation and heating effect on the surface of the electric conductor.[43] The degradation of DWCNTs might be more severe when corona discharge happens because plasma, free radical, as well as ozone can easily destroy DWCNTs at the surface of the discharge electrode. TiO2 is a natural dioxide to wrap the in-layer DWCNTs for avoiding degradation.[28,43] It can overcome the oxidation of DWCNTs efficiently especially during the discharge procedure. DWCNTs and TiO2 composites have shown enhanced thermal conductivity.[50] It appears that the TiO2/CNT-coated electrode exhibits better performance as the heating effect might be perfectly resolved. The highest output power as well as energy conversion efficiency of the TiO2/CNT-coated electrode could be verified and explained excellently. Hence, the DWCNTs and TiO2 nano-composites seem to be an awesome combination to realize satisfactory corona discharge effectiveness for CDIW electrical devices.

4. Conclusion

Dip-coated DWCNTs combined with TiO2 sol have been decorated on the tip of corona needles and have presented scores of merits including lower corona onset voltage, higher CDIW velocity, larger output power, and better energy conversion efficiency. The corona discharge characteristics adhere to the Townsend’s discharge rules as well as the equations presented in the previous publications. In addition to the phenomenological explanation of our experimental data, a series of parameters such as ionic mobility, Reynolds number, and so on for various experimental electrodes are calculated, linearly fitted, and discussed for comparison. The TiO2/CNT-coated device shows the optimal capability for CDIW appliance. The superiority of the coated DWCNTs and TiO2 layers is mainly ascribed to the natural oxidation and protection of TiO2 as well as the particular aspect ratios of DWCNTs. TiO2 and DWCNTs nano-composites may exhibit better electrical and thermal conductivity, less heating effect, and descended contact resistance. In conclusion, DWCNTs with unique high aspect ratios and electrical conductivity companied with natural anti-oxidation of TiO2 contribute to marvelous CDIW manifestation of our apparatus. It is worth drawing on the experience of our neoteric decoration techniques and inspiration for future corona discharge as well as aerodynamics researches.

Acknowledgment

The authors thank Miss. Xiao Meng for her help in figure drawing.

Reference
[1] Robinson M 1962 Am. J. Phys. 30 366 https://doi.org/10.1119/1.1942021
[2] Wu Y Li J Ye J Chen X Li H Huang S Zhao R Ou-Yang W 2017 J. Phys. D: Appl. Phys. 50 395304 https://doi.org/10.1088/1361-6463/aa81a9
[3] Chen I Y Guo M Z Yang K S Wang C C 2013 Int. J. Heat Mass Transf. 57 285 https://doi.org/10.1016/j.ijheatmasstransfer.2012.10.015
[4] Kim B Lee S Lee Y Kang K H 2012 J. Electrostat. 70 438 https://doi.org/10.1016/j.elstat.2012.06.002
[5] He X Jia G Tatsumi E Liu H 2016 Innov. Food Sci. Emerg. 34 135 https://doi.org/10.1016/j.ifset.2016.01.006
[6] Choi S Puligundla P Mok C 2017 LWT — Food Sci. Tech. 75 323 https://doi.org/10.1016/j.lwt.2016.08.063
[7] Wen T Y Wang H C Krichtafovitch I Marnishev A V 2015 J. Electrostat. 73 117 https://doi.org/10.1016/j.elstat.2014.11.002
[8] Yang Z Song H M Jin D Jia M Wang K 2018 Chin. Phys. B 27 085205 https://doi.org/10.1088/1674-1056/27/8/085205
[9] Timmermann E Prehn F Schmidt M Höft H Brandenburg R Kettlitz M 2018 J. Phys. D: Appl. Phys. 51 164003 https://doi.org/10.1088/1361-6463/aab48b
[10] Drew D S Pister K S J 2017 Micromachines 8 141 https://doi.org/10.3390/mi8050141
[11] Audier P Fénot M Bénard N Moreau E 2016 Appl. Phys. Lett. 108 084103 https://doi.org/10.1063/1.4942606
[12] Azooz A A Talal S K 2011 Chin. Phys. Lett. 28 115202 https://doi.org/10.1088/0256-307X/28/11/115202
[13] Eifert A Baier T Hardt S 2013 Appl. Phys. Lett. 103 023114 https://doi.org/10.1063/1.4813556
[14] Lee J R Lau E V 2017 Appl. Therm. Eng. 114 554 https://doi.org/10.1016/j.applthermaleng.2016.12.030
[15] Zhang Y Liu L Chen Y Ouyang J 2015 J. Electrostat. 74 15 https://doi.org/10.1016/j.elstat.2014.12.008
[16] Moreau E Audier P Benard N 2018 J. Electrostat. 93 85 https://doi.org/10.1016/j.elstat.2018.03.009
[17] Chen S vandenBerg R G W Nijdam S 2018 Plasma Sources Sci. T. 27 055021 https://doi.org/10.1088/1361-6595/aabd5f
[18] Li H Jiang L Guo C Zhu J Jiang Y Xiao W Fang C Chen Z 2017 J. Electrostat. 86 59 https://doi.org/10.1016/j.elstat.2017.01.013
[19] Nahan K S Alvarez N Shanov V Vonderheide A 2017 J. Am. Soc. Mass Spectrom 28 2408 https://doi.org/10.1007/s13361-017-1774-0
[20] Bo Z Yu K Lu G Mao S Chen J Fan F 2010 Environ. Sci. Tech. 44 6337 https://doi.org/10.1021/es903917f
[21] Bo Z Yu K Lu G Cui S Mao S Chen J 2011 Energy Environ. Sci. 4 2525 https://doi.org/10.1039/c1ee01140e
[22] Yang W Zhu R Zong X 2016 Nanoscale Res. Lett. 11 1 https://doi.org/10.1186/s11671-015-1209-4
[23] Lv W Q Wang L P Cao Y Z Gu Z W Wang X F Yan Y D Yu F L 2016 Chin. Phys. Lett. 33 095201 https://doi.org/10.1088/0256-307X/33/9/095201
[24] Ye S Liang J Song X Luo S Liang J 2016 Biosyst. Eng. 150 123 https://doi.org/10.1016/j.biosystemseng.2016.07.010
[25] Sadeghian R B Islam M S 2011 Nat. Mater. 10 135 https://doi.org/10.1038/nmat2944
[26] Xu J Xu P Guo P Ou-Yang W Chen Y Feng T Piao X Wang M Sun S 2015 RSC Adv. 5 21755 https://doi.org/10.1039/C4RA16095A
[27] Song Y Li J Ou-Yang W 2019 IEEE Trans. Electron Dev. 66 716 https://doi.org/10.1109/TED.2018.2878738
[28] Xu J Xu P Ou-Yang W Chen X Guo P Li J Piao X Wang M Sun Z 2015 Appl. Phys. Lett. 106 073501 https://doi.org/10.1063/1.4909552
[29] Wu Y Li J Ye J Song Y Chen X Huang S Sun Z Ou-Yang W 2017 J. Alloy. Compd. 726 675 https://doi.org/10.1016/j.jallcom.2017.08.026
[30] Ou-Yang W Weis M Taguchi D Chen X Y Manaka T Iwamoto M 2010 J. Appl. Phys. 107 124506 https://doi.org/10.1063/1.3449078
[31] Dau V T Dinh T X Bui T T Tran C D Phan H T Terebessy T 2016 Exp. Therm. Fluid Sci. 79 52 https://doi.org/10.1016/j.expthermflusci.2016.06.023
[32] Dau V T Tran C D Dinh T X Dang L B Terebessy T Bui T T 2018 IEEE Trans. Dielectr. Electr. Insul. 25 900 https://doi.org/10.1109/TDEI.2018.006881
[33] Gradov O V Gradova M A 2016 Surf. Eng. Appl. Electrochem. 52 117 https://doi.org/10.3103/S1068375516010063
[34] Robinson M 1961 Trans. Am. Inst. Elect. Eng. Part. I: Commun. Electron. 80 143 https://DOI.org/10.1109/TCE.1961.6373091
[35] Bandurin D A Torre I Kumar R K Shalom M B Tomadin A Principi A Auton G H Khestanova E Novoselov K S Grigorieva I V Ponomarenko L A Geim A K Polini M 2016 Science 351 1055 https://doi.org/10.1126/science.aad0201
[36] Liu N Liu J B Yao K L 2017 AIP Adv. 7 125117 https://doi.org/10.1063/1.5001356
[37] Rathi S Lee I Kang M Lim D Lee Y Yamacli S Joh H Kim S Kim S Yun S J Choi S Kim G 2018 Sci. Rep. 8 7144 https://doi.org/10.1038/s41598-018-22355-0
[38] Ou-Yang W Mitoma N Kizu T Gao X Lin M F Nabatame T Tsukagoshi K 2014 Appl. Phys. Lett. 105 163503 https://doi.org/10.1063/1.4898815
[39] Nakayama K Ou-Yang W Uno M Osaka I Takimiya K Takeya J 2013 Org. Electron. 14 2908 https://doi.org/10.1016/j.orgel.2013.08.002
[40] Ladouceur H D Baronavski A P Lohrmann D Grounds P W Girardi P G 2001 Opt. Commun. 189 107 https://doi.org/10.1016/S0030-4018(01)01012-4
[41] Housley D Huddle T Lester E Poliakoff M 2016 Chem. Eng. J. 287 350 https://doi.org/10.1016/j.cej.2015.11.013
[42] Zhang J Deng Y Hao T T 2018 Chin. Phys. B 27 128101 https://doi.org/10.1088/1674-1056/27/12/128101
[43] Zhu W D Wang C W Chen J B Li D S Zhou F Zhang H L 2012 Nanotech. 23 455204 https://doi.org/10.1088/0957-4484/23/45/455204
[44] Chen X Jiang T Sun Z Ou-Yang W 2015 Appl. Phys. Lett. 107 114103 https://doi.org/10.1063/1.4931463
[45] Shin D H Baek S H Ko H S 2016 Int. J. Heat Mass Transf. 93 516 https://doi.org/10.1016/j.ijheatmasstransfer.2015.10.029
[46] Zhou D Tang J Kang P Wei L Zhang C 2018 J. Electrostat. 96 99 https://doi.org/10.1016/j.elstat.2018.10.007
[47] Alivov Y Klopfer M Molloi S 2010 Appl. Phys. Lett. 96 243502 https://doi.org/10.1063/1.3454244
[48] Zhang G Huang S Y 2013 Physics 42 100 https://doi.org/10.7693/wl20130203
[49] Zhou K K Xu N Xie G F 2018 Chin. Phys. B 27 026501 https://doi.org/10.1088/1674-1056/27/2/026501
[50] Khosravifard E Salavati-Niasari M Dadkhah M Sodeifian G 2012 J. Nanostruct. 2 191 https://DOI.org/10.7508/JNS.2012.02.006